The nonrelativistic limit of Dirac-Fock codes: the role of Brillouin configurations 
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We solve a long standing problem with relativistic calculations done with the widely used Multi- 
Configuration Dirac-Fock Method (MCDF). We show, using Relativistic Many-Body Perturba- 
tion Theory (RMBPT), how even for relatively high-Z, relaxation or correlation causes the non- 
relativistic limit of states of different total angular momentum but identical orbital angular momen- 
tum to have different energies. We show that only large scale calculations that include all single 
excitations, even those obeying the Brillouin's theorem have the correct limit. We reproduce very 
accurately recent high-precision measurements in F-like Ar, and turn then into precise test of QED. 
We obtain the correct non-relativistic limit not only for fine structure but also for level energies and 
show that RMBPT calculations are not immune to this problem. 

PACS numbers: 31.30.Jv, 32.10.Fn, 31.25.Eb 
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Relativistic atomic structure codes, mostly MCDF 
packages, are now of widespread use in many sectors 
of physics, and the need for reliable, relativistic calcu- 
lations is stronger than ever (see, e.g., Q for examples 
in Astrophysics). However, the difficulties of doing reli- 
able calculations are numerous, and still largely under- 
estimated. For example a puzzle that was noted already 
twenty-two years ago [2| has never been solved, although 
it may lead even in very simple calculations to wrong en- 
ergy values. In Ref. Q it was shown that relativistic 
self-consistent field procedures do not produce, in a num- 
ber of cases, the correct non-relativistic limit of zero for 
the fine structure splitting (FSS) when the speed of light 
is tuned to infinity. Ref. |2| suggested as a remedy ex- 
plicit calculation of this non-relativistic offset (N.R.) and 
subsequent subtraction of it from the relativistic result, 
although no justification for the procedure was provided. 
Moreover this paper said nothing on how to correct indi- 
vidual energy levels. Here we will penetrate the origin of 
the non-relativistic shift using the tools of perturbation 
theory and advanced MCDF calculations. We use these 
tools to show the role of relaxation in the N.R. offset, 
and prove that the inclusion of specific mono-excitations 
in the MCDF basis removes it. We also provide justifica- 
tion to the subtraction procedure and show that not only 
the FSS need to be corrected, but also the level energy, 
e.g., when transitions between different shells are stud- 
ied. It is also worth noting that this problem appears in 
the Optimized Level (OL) scheme when each level energy 
and wavefunction is optimized separately. This scheme 
is used only when the highest accuracy for correlation is 
required. Often the average (AL) level scheme is used, in 
which the same J-average wavefunction is used to calcu- 
late the energy of all FS component. In the AL scheme 
the N.R. offset does not appear, but the accuracy is much 
lower. 



We will concentrate on the ground state configuration 
of a F-like ion which was used as a model system already 
in Ref. Q as accurate measurements have been performed 
very recently Q . With high experimental accuracy, even 
for Z = 18, it is important to be aware of this problem 
which seriously affects the comparison with experiment 
on the present day level. We will further present accurate 
calculations of the fine structure splitting in F-like argon 
both with Relativistic Many-Body Perturbation Theory 
(RMBPT) and with the Multi-Configuration Dirac-Fock 
(MCDF) method. It is shown that by comparison with 
accurate experimental results |3| it is possible to test the 
calculations on self-energy and other radiative corrections 
in a true many-electron surrounding. 

With RMBPT the fine-structure splitting in a F-like 
system is calculated as the binding energy difference be- 
tween the 2pi/2 and 2p 3 / 2 electron in the correspond- 
ing Ne-like system. The lowest order approximation of 
this binding energy is the negative of the orbital energy 
of the removed electron in the Hartree-Fock approxima- 
tion. The remaining electrons are at this stage considered 
as frozen in their orbitals in spite of the removal of one 
electron. The most important correction to this first ap- 
proximation is the relaxation of the electrons due to the 
presence of the hole. The term relaxation usually denotes 
the correction found by a single configuration restricted 
Hartree-Fock (or Dirac-Fock in the relativistic case) cal- 
culation in the presence of the hole. The non-relativistic 
shift has its origin already at this level and we will now 
concentrate on this shift and postpone the discussion of 
higher order corrections. 

To analyze the relaxation for a one-hole state with per- 
turbation theory it is natural to start from the closed 
shell system and systematically correct for the removal 
of one electron. Fig. ^ shows the contributions entering 
in second order. Fig. ^a-b) show fluctuations to two 



2 



ht 



c T Ts a 



ht h t 

(a) (b) 




decoupling of spin and orbital angular momentum can- 
not be done without summing over all total angular mo- 
menta, j. An unambiguous way to see how this influences 
our example of Fig.^a) with orbital a = h and i s = it, is 
to compare the angular contribution non-relativistically 
and relativistically. The electron-electron interaction is 
expressed as 



FIG. 1: Illustration of the second order energy contributions 
to a one-hole system. Diagrams(a-b) show fluctuations to two 
holes and one excited orbital and diagrams (c-d) double ex- 
citations (correlation). Downgoing single arrows denote core 
orbitals, downgoing double arrows denote the hole and upgo- 
ing arrows denote excited orbitals. 



holes and one excited orbital and Fig. ^c-d) true double 
excitations. The relaxation, i.e., the effects included by 
a single configuration restricted Hartree-Fock calculation 
is in perturbation theory part of Fig. ^a-b); the ones 
where the hole is not fluctuating and the excitation from 
an orbital preserves its angular symmetry. The lowest or- 
der relaxation correction to an orbital b can consequently 
be written 
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where h denotes the removed electron, the curly brack- 
ets antisymmctrization, V\2 the two-electron interaction, 
and the minus sign is due to the removal of h. The energy 
corrections are then calculated as 
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In this way all types of diagrams in Fig. ^ (a-b) with ei- 
ther orbital a = h (and b=c) or c = h (and b=a) and 
i s = ib, js = jb are included, i.e., the single excitations 
that preserve the angular structure. It can be noted that 
these single excitation contributions form a class of di- 
agrams that can be summed until convergence in an it- 
erative scheme, see, e.g., Refs. 0,0. Here we will not 
pursue this line, however, since our purpose is to ana- 
lyze the relaxation in the non-relativistic limit and show 
why a state with a hole in nij = i-\/2 and one with a hole 
nij— (.+1/2 do not reach the same energy in this limit. For 
this it is sufficient to study relaxation in second order. 

As an example, take diagram Fig. Q[a) with orbital 
a = h and i s — ib, a typical relaxation contribution. The 
orbitals used to evaluate the diagram are solved using 
the Hartree-Fock potential from the closed shell core and 
the radial part of the 2p^/ 2 and the 2p 3 / 2 orbital will be 



identical when we let c — > oo. The problem comes instead 
from the spin- angular part. Since 
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where k denotes the rank of the spherical tensor operator 
C, which works on the orbital part of wavefunctions. Non 
relativistically the angular part can be evaluated as 
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This is in fact identical to the following expression in the 
coupled space where two extra sums appear over inter- 
mediate total angular momenta 
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That these two expressions give the same result can 
be understood by standard angular momentum algebra 
techniques. 

In a restricted Dirac-Fock calculation there will be no 
sums over intermediate angular momenta. Instead only 
jh = jh' is allowed, i.e., the hole is not allowed to fluctu- 
ate to the other fine structure component, and jV = % is 
required , i.e., the corrections to orbital 6 do not change 
its angular structure. The spin-angular part used is thus 
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which will clearly not produce the same result as Eq. J5j, 
and which further cannot give identical results for, e.g., 
jh = 1/2 and jh — 3/2, which is easily seen from the 
k = 2 contribution which is zero for jh — 1/2, but not 
for jh — 3/2. The difference can also be readily demon- 
strated numerically for a system as F-like neon where 
the second order contribution, Eq. [21 to the relaxation 
gives an unphysical fine structure offset of 0.024 eV in 
the c — > oo limit. Following the recipe from Ref. 0] and 
correcting the result calculated with the true value of c 
with this offset, we obtain a relaxation contribution to 
the fine structure splitting of —0.058 eV. After iteration 
of the relaxation contributions 0, the corrected value 
reaches ~ —0.050 eV, in line with the MCDF Coulomb re- 
laxation contribution of ~ —0.049 eV, listed in the third 
section of Table [I] This value has been corrected using 
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TABLE I: Summary of the contributions. All calculations 
use the 2002 values for fundamental constants [(J (eV). 
Experimental values are from wavelength provided in Ref. 0] 
converted to vacuum values using y| 





2P1/2 


2P3/2 


A 


Contributions 






Ne-like DF orb. ener. 


426.50002 


424.13211 


2.36791 


A DF-Breit 


-0.22659 


-0.13576 


-0.09083 


h.o. retardation 


-0.00011 


0.00079 


-0.00090 


QED corr. 


0.01353 


0.00755 


0.00598 


Contributions specific to RMBPT 




2nd order core-core, Coul 


-4.48509 


-4.42587 


-0.05921 


core-core, Breit 


-0.01187 


-0.00814 


-0.00373 


correlation, Coul 


2.56726 


2.55763 


0.00962 


correlation, Breit 


0.02391 


0.02018 


0.00373 


h.o. contr. (Coul.+Breit) 


0.16559 


0.15885 


0.00674 


A DF Breit orbitals 


0.00198 


0.00043 


0.00156 


Total (RMBPT) 


424.54863 


422.30777 


2.24086 


Experiment 






2.24010 


Contributions specific to MCDF (N.R. offset subtracted) 


Relaxation (Coul) 


-3.10800 


-3.05931 


-0.04869 


Relaxation (Breit) 


-0.00406 


-0.00314 


-0.00092 


Correlation (Coul — > 5g) 


1.42466 


1.39604 


0.02862 


Correlation (Breit — ► 5g) 


-0.01359 


0.00741 


-0.02100 


Total (MCDF) 


424.58585 


422.34569 


2.24016 


Experiment 






2.24010 



the same procedure. The small difference is probably due 
to small differences in the classification of relaxation and 
correlation contributions. The lesson here is that since 
the summation over all possible couplings of spin and 
orbital angular momenta of the intermediate states are 
necessary to reproduce the uncoupled situation a correct 
non-relativistic limit cannot be achieved with any single 
configuration self consistent field calculation. In other 
words, still for the system under consideration, one has 
to include more than one configuration relativistically to 
reproduce the single configuration non-relativistic result 
in a relativistic framework. In RMBPT the full contri- 
bution from Fig. ^a-b) produces no N.R. offset, but any 
attempt to speed up the convergence of the perturbation 
expansion by singling out the important subclass that 
only involve energy denominators as in Eq. ^ will do so. 
The RMBPT results shown in Tabled are obtained with- 
out any such procedure and has a correct non-relativistic 
limit by construction. 

With several configurations included it should in prin- 
ciple be possible to reach the correct non-relativistic 
limit, in practice one can, however, generally not achieve 
this in a truncated calculation. In practice the number of 
configurations has to be truncated for all but the smallest 
systems. It is common to truncate after double excita- 



TABLE II: Contributions to the MCDF energy affected by 
the N.R. offset (eV). "AE doub. Exc. -tn = i": correlation 
energy for the configuration space which include all double 
excitations up to principal quantum number n — i. Rel. Val.: 



Relativistic Value. N.R. 


Off.: Offset 


obtained at 


the non- 


relativistic limit. 










Rel. Val. 


N.R. Off. 


Diff. 


Dirac-Fock Coulomb 


2.31626 


-0.00148 


2.31774 


Brillouin sing 


;le excitations excluded 




AE Exc. -> n = 3 


-0.01855 


-0.02086 


0.00231 


AE Exc. -> n = 4 


-0.01421 


-0.01926 


0.00505 


AE Exc. -> n = 5 


-0.01641 


-0.02247 


0.00606 


Total 


2.21621 


-0.02395 


2.24016 


Diff. With Exp. 


-0.02389 




0.00006 


All single and double excitations included 




AE Exc. -> n = 3 


-0.00371 


-0.00582 


0.00211 


AE Exc. — > n = 4 


0.00445 


-0.00037 


0.00482 


AE Exc. -> n = 5 


0.00661 


0.00075 


0.00586 


Total (S.E.S. Welton) 


2.23923 


-0.00073 


2.23996 


Diff. With Exp. 


-0.00087 




-0.00014 



tions from the dominating configuration, but just as dou- 
ble excitations are needed to be added to the single exci- 
tations to obtain the correct non-relativistic limit, triple 
excitations will be needed to be added to correspond- 
ing double excitations and so on. Since higher multiple 
excitations are less important the remaining offset will 
however decrease steadily. 

We now proceed to demonstrate the vanishing of the 
non-relativistis offset in an essentially complete MCDF 
calculation. In the present calculation we have added to 
the original configuration all single and double excita- 
tions up to a given maximum n and I. Note that one 
has to be careful in considering the meaning of single 
and double excitations. For example the ls 2 2s 2 2p 4 3p is 
a single excitation in the LS coupling sense. Yet in jj 
coupling it gives rise to 5 configurations in the J = 1/2 
case, two of which are double excitations in the jj sense 
(2pi/ 2 2p3 /2 -> 2p 2 1/2 2pj /2 3p 1/2 and 2p 2 1/2 2p 2 3/2 3p 3/2 ).We 
went from 3c? to 5g for the case with a normal speed of 
light, and up to 6h for the non-relativistic limit. This 
represents respectively 299, 1569, 4339 and 9127 fully re- 
laxed jj configurations for the J — 1/2 case, and 456, 
2541, 7356 and 15915 for the J = 3/2. The calculations 
are repeated with different lists of configurations. In one 
group of calculations, we include all single and double 
excitation in the jj sense, except for the "Brillouin sin- 
gle excitations", i.e., those that should contribute only in 
third order, as stated by Brillouin's theorem [9. fiol Hlj. 
These excitations are often excluded since they compli- 
cate the numerical convergence. Again we use here Bril- 
louin's theorem in the jj sense, i.e., we exclude all con- 
figurations transformed from the initial one by replacing 
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FIG. 2: Comparison of the non-relativistic offset for Be-like 
and F-like argon, evaluated including all single excitations, or 
only those not obeying Brillouin's theorem. 



space is used, but only if all single configurations are in- 
cluded. In practice excluding Brillouin single excitations 
and then subtracting the N.R. offset leads to the same 
value, but numerical convergence of the self-consistent 
field process is much easier in the latter case. Finally, 
failing to account for the N.R. offset leads to poor re- 
sults, even at a moderately large Z , a fact that may not 
have received enough attention in many MCDF calcula- 
tions. The present work also shows that similar prob- 
lems can happen in RMBPT calculations if subclasses of 
important effects are singled out and by themselves are 
treated to higher order. The improved convergence will 
then come at the expense of an N.R. offset. This fact had 
not been recognized before. 



an orbital with quantum numbers n, n by one with n', k, 
where k is the Dirac angular number. In a second group 
we include all single and double excitations. In both 
groups, we do calculations once with only the Coulomb 
interaction between electrons used in the evaluation of 
wavefunctions and energies, and once with the full Breit 
interaction in the evaluation of wavefunctions and mix- 
ing coefficients. This allows to include high-orders of the 
Breit interaction in the calculation. In each group the 
Coulomb only calculation is done also a second time with 
a large value for the speed of light. The evolution of the 
N.R. shift as a function of the maximum excitation used 
in the MCDF process is plotted in Fig. |3 for both F- 
like and Be-like ions, to show the generality of what is 
observed: the N.R. offset tends to a non-zero constant 
value when Brillouin configurations are excluded, and to 
zero when all single excitations are included. 

The different contributions to the MCDF result, 
and the variation of the correlation energy and non- 
relativistic offset with and without Brillouin configura- 
tions are presented in Table ITU When comparing both 
results, it is clear that excluding Brillouin single exci- 
tations and then subtracting the N.R. offset leads to the 
same result as including the Brillouin configurations. The 
agreement with experiment and with RMBPT results is 
excellent in both cases, even though the quality of the 
convergence when including all single excitations is not 
as good as when Brillouin ones are excluded. Moreover, 
the inclusion of all single excitations enables also to cor- 
rect the energy of a level as shown in Table IIIII which 
was not possible with the subtraction method. Finally 
we note that the evaluation of the radiative corrections, 
the self-energy screening (SES) with the help of the Wel- 
ton approximation [12[ leads to a very good agreement 
with experiment. 

In conclusion we have proven, by comparing RMBPT 
and MCDF results, that the N.R. offset is due to relax- 
ation and should go away when doing a complete calcu- 
lation. We then showed that in the MCDF case, the 
offset is going to zero if a large enough configuration 

TABLE III: Change in the ground state (J = 3/2) correlation 
energy due to Brillouin single excitation 

Conf. No-brillouin All single Diff 

corr -» 3d -5.1792 -5.1989 -0.0196 

corr -> 4/ -7.7349 -7.7603 -0.0255 

nnrr _>. Kr, _S R ^ 5 1 _S RS71 _fl fl^n 
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